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The surface density of the star formation rate in different galaxies, as well as in 
different parts of a single galaxy, scales nonlinearly with the surface density of the total 
gas. This observationally established relation is known as the Kennicutt-Schmidt star 
formation law. The slope of the star formation law has been shown to change with the 
density of the gas against which the star formation rate is plotted. This dependence 
implies a nonlinear scaling between the dense gas and the total gas surface densities 
within galaxies. Here, we explore a possible interpretation of this scaling as a property 
of the geometry of the interstellar medium (ISM), and we find that it arises naturally if 
the topology of the ISM is multifractal. Under the additional assumption that, at very 
high densities, the star formation timescale is roughly constant, the star formation law 
itself can also be recovered as a consequence of the multifractal geometry of the ISM. 
The slope of the scaling depends on the width of the global probability density function 
(PDF), and is between 1.5 and 1.6 for wide PDFs relevant to high-mass systems, while 
it is higher for narrower PDFs appropriate for lower-mass dwarf galaxies, in agreement 
with observations. 
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1 INTRODUCTION 

The dependence of the star formation rate on gas density 
on large scales is a subject of intense observational and the- 
oretical investigation, as it is both an important input for 
models of galaxy formation and evolution and chemical evo- 
lution calculations, as well as a critical test for theories of 
interstellar medium evolution and star formation. Let us de- 
fine the surface density of the star formation rate, Ssf, as 
the mass of gas being converted into stars per unit time per 
unit surface area, and E gas as the total gas mass per unit 
surface area. The two quantities have been found to obey a 
non-linear, power-law relation: 

E S f oc Egir , (1) 

where n gas w 1.5 for a large range of surface densities 
and system morphologies. This correlation is known as the 
Kennicutt-Schmidt law of star formation (Schmidt 1959, 
Kennicutt 1989). The star formation law has been estab- 
lished through observations of the global star formation and 
gas density in different galaxies (e.g. Kennicutt 1998, Misiri- 
otis et al. 2004, Komugi et al. 2005); of local gas and star 
formation densities in different galaxies (e.g. Wong and Blitz 
2002, Boissier et al 2003); and of local gas and star forma- 
tion densities within a single galaxy (e.g. Misiriotis et al. 



2006 in the case of the Milky Way, Schuster et al. 2007 in 
the case of M51). 

Traditionally, this correlation has been interpreted in 
the literature as a result of the density dependence of the 
star-formation timescale. The star formation rate density 
can be expressed as a ratio of the density of the gas available 
for star formation over a timescale relevant to the conver- 
sion of the available gas to stars. Examples of such timescales 
that have been suggested in this framework are the free-fall 
timescale, the turbulence crossing time over the scale height 
of the galaxy, the collapse timescale of large expanding shells 
with low Mach number (see e.g. Elmegreen 2002b and ref- 
erences therein), the orbital timescale of the galactic disk 
(Silk 1997), and the timescale for gas accumulation along 
magnetic flux tubes parallel to spiral arms, into the valleys 
created by the magnetic Rayleigh- Taylor (Parker) instabil- 
ity (see e.g. Shu et al. 2007). These timescales all scale as the 
inverse square root of density, resulting in an overall scaling 
of the star formation density as gas density to the 1.5. A 
more elaborate treatment based on the same principle was 
presented by Krumholz and McKee (2005). This interpre- 
tation of the star formation law is very tempting because 
it is conceptually simple and elegant, and because it gives 
the same result for several different processes which may be 



2 K. Tassis 



controlling the star formation timescale. However, a series 
of more recent observations have indicated that this simple 
picture may not be a complete interpretation of the star 
formation law. 

First of all, the value n gas ~ 1.5, although consistent 
with observations of a large range of star-forming systems, 
is by no means unique. Appreciable scatter exists that can- 
not be entirely attributed to observational uncertainties and 
cannot be comfortably reconciled with a theoretically set 
single value of n gas . For example, if only gas with local vol- 
ume densities representative of molecular cloud cores is in- 
cluded in the gas surface density (as is the case when tracers 
sensitive to higher density gas are used, such as HCN) the 
resulting scaling is closer to linear, 

Ssf cx , (2) 

where Edcn, gas is the surface density of the denser gas, and 
^dense is now closer to 1 rather than 1.5 (e.g. Gao and 
Solomon 2004, Wu et al. 2005). Observations of dwarf galax- 
ies on the other hand seem to favor much steeper scalings 
(values of n gas much larger than 1.5). For example, de Blok 
and Walter (2006) find for NGC 6822 a slope 2.04; Heyer et 
al. (2004) find for M33 a much steeper slope, equal to 3.3; 
in the case of IC10, Leroy et al. (2006), find a slope of 1.4 
for molecular gas and a slope much steeper than 1.5 for the 
total gas density (see their Fig. 14). 

In addition, since star formation is an inherently local 
phenomenon, the connection between the small scales (~ 
0.1 pc) where stars form, and the large scales (~ 1 kpc) that 
the gas density observations and the associated timescales 
refer to, may be quite complicated. The traditional inter- 
pretation implicitly assumes that a meaningful mean density 
can be assigned to the large scales where £ gas is determined. 
However, this need not be the case. As a counter-example, 
we can consider the situation where several dense objects are 
sparsely distributed over a region of ~ kpc dimensions. The 
mean density we would derive would be rather low, but it 
would not be representative of any object in the region. Con- 
sequently, the star-formation timescale calculated at that 
mean density would not represent the real timescale over 
which any of the objects in the region evolve. Instead, the 
relevant timescale for star formation should be constructed 
out of the densities of the individual star- forming sites rather 
than the large-scale mean density. 

Finally, the interstellar medium (ISM) is a highly com- 
plex system, and its structure and properties are the result 
of a synthesis of nonlinear phenomena and instabilities, not 
necessarily connected with one another. Star formation on 
the other hand only occurs at the density peaks of the ISM, 
and it is likely that the star formation law may be con- 
veying information about the connection between the large- 
scale low-density ISM and its small-scale high-density peaks, 
rather than information about the star formation process it- 
self. 

The connection between the high-density and low- 
density ISM gas is imprinted in different incarnations of the 
star formation law. Observations have established that the 
slope of the scaling of the star formation rate surface density 
with the gas surface density depends on the minimum den- 
sity of the gas accounted for in observations [Eqs. |TJ and 
{2}]. If now we demand that Eqs. ([!]) and Q both hold, this 



implies a scaling between the surface densities of total gas 
and dense gas, 

^den, ga s CX S gas , (3) 

with n = n gas /ndcnso ~ n gas w 1.5. 

In this paper we investigate this connection, and we 
examine the extent to which such a nonlinear scaling of the 
dense gas surface density with the total gas surface density 
can be interpreted as a topological property of the ISM. In 
particular, we investigate whether such a scaling can arise 
naturally in multifractal topologies. 

While fractal structures are characterized by a unique 
non-integer generalized dimension, a spectrum of general- 
ized fractal dimensions is required to fully describe multi- 
fractal structures. In a multifractal structure, we can define 
a scaling exponent which characterizes the structure locally, 
around each point. Sets of points that share the same scal- 
ing exponent constitute a fractal. The distribution of fractal 
dimensions of all such fractals is the multifractal spectrum, 
which describes the distribution of geometries present in a 
complex structure. 

The gas and dust in the ISM is organized in compli- 
cated structures with irregularities present over a wide range 
of scales. Different studies have used different techniques to 
quantify the properties of these structures and describe the 
geometry of the ISM: e.g., autocorrelation function (Dick- 
man and Kleiner 1985), power spectra (Stutzki et al. 1998); 
perimeter/area relation (Bazell and Desert 1988; Scalo 1990; 
Vogelaar and Wakker 1994; Westpfahl et al. 1999 in HI for 
members of the M81 group), wavelet analysis (Langer et 
al. 1993), principal component analysis (Heyer and Schlo- 
erb 1997), multifractal scaling (Chappell and Scalo 2001). 
These studies find that structure in the ISM has a scale- 
free, hierarchical appearance. This complexity over a wide 
range of scales has lead to the proposition that a multifractal 
ISM geometry with a lognormal probability density function 
(PDF) is consistent with many of the observed features, in- 
cluding the mass and size distribution of clouds and the 
ISM filling factor (Elmegreen 2002a; Elmegreen et al. 2006). 
The physics is usually associated with compressible turbu- 
lence (e.g., Mac Low and Klessen 2004, Scalo and Elmegreen 
2004), however the multifractal scaling is much more generic 
and may arise from many different multiplicative hierarchi- 
cal processes (e.g. Colombi, Bouchet, and Schaeffer 1992; 
Sylos Labini and Pietronero 1996). 

It should be noted that the ISM cannot be strictly scale- 
free, since optical observations of a variety of galaxies sug- 
gest an upper characteristic scale (~1 kpc) that can either 
be the Jeans length (e.g. Elmegreen 2004) or the wavelength 
of the Parker instability (Mouschovias, Shu and Woodward 
1974). A lower characteristic scale also exists, at the or- 
der of ~0.1 pc (Blitz and Williams 1997; Barranco and 
Goodman 1998; Goodman et al. 1998). However, in between 
these scales, a variety of different processes and instabilities 
can lead to multifractal scalings and lognormal PDFs (e.g. 
Vazquez- Semadeni 1994; Kravtsov 2003; Wada & Norman 
2007). 

The purpose of this work is to examine whether a non- 
linear scaling of the dense gas surface density with the to- 
tal gas surface density similar to that of Eq. © can arise 
naturally in multifractal ISM geometries, whether it is ro- 
bust enough to resemble the persistence and uniformity of 
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rigaa/Hdense over a variety of star forming systems with dif- 
ferent local conditions, and whether it can allow for different 
slopes seen e.g. in observations of dwarf galaxies under phys- 
ical conditions similar to those found in these systems. 

A geometrical explanation to Eq. ((3j has an additional 
tantalizing potential consequence. Equation ((2)1 implies that 
the star formation rate surface density scales linearly with 
the very dense gas, or, equivalently, that at those high den- 
sities the star formation timescale is roughly constant and 
independent of Sden.gas: 

E sp = Sdcn.g- ( 4) 
T"SF, dense 

with tsf, dense ~ const. However, in this case Esf oc Eden.gas 
and if E gas oc Sj on gas due to geometry, then Es F oc Eg as . 
The observationally motivated assumption of a star forma- 
tion timescale independent of Eden at high densities, in com- 
bination with a topologically-driven scaling between dense 
and total gas, can thus provide a geometrical interpretation 
of the star formation law itself. 

This paper is structured as follows. Our formulation and 
our algorithms for the reproduction of multifractal three- 
dimensional geometries are discussed in Section [2] Our re- 
sults are presented in Section [3] We discuss our findings in 
Section [H and we summarize our conclusions in Section \5\ 



2 FORMULATION 

In order to investigate how the mean surface density of very 
dense gas correlates with the surface density of the total gas 
in different multifractal ISM topologies, we need to generate 
a multifractal structure, and perform appropriate "observa- 
tions" of its total and dense "gas" surface densities. In this 
section, we describe these procedures. 

To generate a multi-fractal, we adapt the algorithm of 
Borgani et al. (1993), which is a modification of the j3 model 
(Frisch et al. 1978) and the random-/3 model (Benzi et al. 
1984) of fully developed turbulence. We start with a parent 
three-dimensional cube of side L, which we divide into 2 3 
equal- volume subcubes, each of which inherits some fraction 
fi of the parent-cube mass, where i = 1, .., 8. Mass conserva- 
tion implies that J^? =1 fi = 1- We repeat the fragmentation 
(where each subcube now becomes a parent-cube) H times. 
The fraction of the total mass contained in each final cube 
of volume L 3 /2 3H depends on its fragmentation history. 

The properties of the final structure depend on the 
choice of fi. If all fi are non-zero and equal, then the re- 
sult is a homogeneous structure. If some of the fi are zero 
and all of the non-vanishing fi are equal, then the result is 
a monofractal structure, with a fractal dimension uniquely 
set by the choice of fi. If the non- vanishing /, are not equal, 
then a multi-fractal structure results. To obtain a pattern- 
free structure that still retains the scaling properties of a 
fractal (self-similarity) or a multifractal (self-affinity), the 
pattern of fi is not kept constant in all iterations: each fi 
value is assigned to a random subcube i in each iteration. 

Each multifractal cube is taken to represent one ISM 
topology. We investigate how the total gas surface density 
correlates with the dense gas surface density in different 
parts of this object in the following way: once a multifractal 



Table 1. Mass fractions fi and levels of refinement H for the 
seven multifractals investigated in our study. Multifractals are 
listed here in order of decreasing PDF width (see Fig. |3j. 
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cube is constructed according to the algorithm above, a line- 
of-sight is selected, parallel to one side of the cube. The face 
of the cube which is perpendicular to that line is then split 
into 2 2( - H ~ 3 ^ square patches (each side of the face is split into 
2 H ~ 3 segments). For each of these patches, an "observation" 
is made of its total surface density, by summing up all of the 
mass within the volume defined by the patch and extending 
along the line of sight, and dividing by the surface area of 
the patch. Similarly, the "dense gas" surface density is calcu- 
lated by summing up all of the mass within regions of local 
volume density above the "dense gas" threshold within the 
same volume, and dividing by the surface area of the patch. 

To test the robustness of our results against variations 
of the detailed properties of the multifractal cube, we have 
repeated our measurements for seven distinct multifractals. 
The properties of these multifractals (fi and number of re- 
finement levels, H) are given in Table [1] 



3 RESULTS 

The first questions we seek to answer are whether the scaling 
between the very dense gas and total gas in a multifractal 
medium resembles that of Eq. ((3j); and whether such a scal- 
ing is robust enough against variations in the detailed prop- 
erties of the multifractal (parameterized in our model by the 
values of the fi and H) so that such a scaling can be viewed 
as an intrinsic property of multifractal geometries, rather 
than as a result of fine-tuning of the multifractal model. 

We address these questions in Fig. [T] which shows the 
surface density of dense gas, Eden.gas, in this case defined 
as gas with local volume density at least 3 x 10 3 times the 
mean volume density, as a function of the surface density of 
total gas Egas, for three different multifractals (A, B, and 
C). Each datapoint corresponds to a different patch within 
the multifractal cube. The three fractals depicted in this 
figure have different mass fractions, fi, however they all have 
"wide" global PDFs (see upper panel of Fig.[3]and discussion 
below). In all cases, the scaling between E gas and Eden.gas 
is clearly nonlinear (significantly deviates from the linear 
relation plotted in Fig. [1] with the solid line), and has a 
scaling slope of « 1.5. 

To understand the origin of the nonlinearity of the scal- 
ing, we need to examine how the gas is distributed over dif- 
ferent local densities. This distribution is quantified by the 
PDF, which is defined as the distribution of volume frac- 
tion with respect to local density (where the local density 
is measured in units of the mean density of the cube). In 
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Figure 1. Scaling of dense gas surface density vs. total gas surface 
density for multifractals A (green), B (red), and C (purple). The 
solid and dashed lines correspond to the £den,i 
Eden, gas °c ^gks scalings respectively. 



order to have a nonlinear scaling between the surface densi- 
ties of total gas and gas density peaks, the PDFs of different 
patches within the cube must be different from each other. 
This requirement is straight forward to prove. Let the local 
PDF within a patch i of volume Vi be g%{p) = (l/Vi)dV/dp. 
Then, the total dense gas mass in patch i is an integral of 
the PDF above some high density threshold pden, 



M de 



Vi 



9i{p)pdp, 



(5) 



and the total gas mass of the same patch is an integral of 
the PDF over all densities, 



-A^gas,i — Vi 



9i(p)pdp, 



while the ratio of surface densities is given by 



■+l Ai Jo°° 9i{p)pdp 



(6) 



(7) 



where A; is the surface area of the patch perpendicular to 
the line of sight. If now the PDF remains the same among 
different patches (all gi are identical), the right-hand side of 
Eq. (0 is always constant, resulting in an identically linear 
scaling, E den E as oc Sgas- A similar result is obtained, as long 
as the PDF spread is sufficiently large, if the gi have the 
same spread but different mean densities. A nonlinear scal- 
ing can only arise if different patches exhibit different PDF 
spreads. This feature is an inherent property of multifrac- 
tals, as we show below. 

Figure [2] shows local PDFs for four different patches of 
varying mean density for multifractal A. In our multifrac- 
tal geometry, the local PDFs are not identical in spread to 
each other, and they are not identical to the global PDF 
of the cube (plotted with the thin solid black line in the 
upper panel of Fig. [3|. Instead, they are truncated at dif- 
ferent densities depending on their total gas surface density. 
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Figure 2. Probability density functions of different patches of 
multifractal A. The total gas surface density is lowest for the 
patch represented by the solid black line, and progressively in- 
creases for the patches depicted in dashed red, long-dashed green, 
and dotted blue respectively. 



The total gas surface density is lowest for the patch repre- 
sented by the solid black line, and progressively increases for 
the patches depicted in dashed red, long-dashed green, and 
dotted blue respectively. The lower the total surface den- 
sity of the patch, the lower the density at which the PDF is 
truncated. It is this effect which is the origin of the nonlin- 



earity in the scaling of Eden, gas with E gas . A similar depen- 
dence of the local PDF on the mean density was pointed out 
by Kravtsov (2003) in the case of cosmological simulations, 
where the star formation law was reproduced in large scales, 
despite the fact that the local star formation recipe featured 
a timescale constant with density. 

We have demonstrated that a nonlinear scaling between 
very dense and total gas surface densities arises naturally in 
multifractal geometries, and we have traced its origin in the 
property of the local multifractal PDF to extend to higher 
local densities in regions of higher mean density. However, 
there are still several points to be addressed before we can 
draw conclusions on how robust this scaling really is, and 
how it may relate to the star formation law. These points 
are the following. 

(i) The definition of "dense gas" we used in our discussion 
of Fig. [T] was rather arbitrary. We would like to understand 
how the slope of the scaling may depend on this definition. 

(ii) Although we have shown that different multifrac- 
tals show similar behaviour, we would like to understand 
whether scalings with slope other than « 1.5 can also be 
recovered, and under what conditions. 

(iii) Although the detailed values of ft , together with the 
number of refinement levels, uniquely determine the multi- 
fractal scaling properties, it would be desirable to identify a 
global property of the multifractal, which controls the slope 
of the scaling. 

(iv) Finally, we would like to examine whether the condi- 
tions under which scalings with slopes different that « 1.5 
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Figure 3. Upper panel: probability density functions for each 
multifractal. Different line styles and colours correspond to dif- 
ferent multifractals as follows: thin solid black — ► A; dotted red 
— » B; long-dashed purple — > C; thick solid green — ► D; dashed 
blue — > E; double-dot-dashed cyan — > F; dot-dashed magenta — » 
G. Lower panel: dependence of the exponent n of the scaling be- 
tween dense gas surface density and total gas surface density on 
the value of "dense gas" density threshold (in units of the mean 
density of the multifractal); line styles and colours as in the upper 
panel. 



appear in our simulated multifractals possibly mirror situa- 
tions in nature where the n gas slope of the star formation law 
entering Eq. <(3j also has different values (as, for example, 
in the case of dwarf galaxies) . 

We address these issues with the help of Figs. [3] and [4] 

The upper panel of Fig. [3] shows the PDF of the cube 
as a whole for each multifractal described in Table [T] The 
global PDF can be fitted well by a lognormal in all cases. 
The lower panel of the same figure shows the dependence of 
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Figure 4. Dependence of scaling slope n on the definition of 



S gas . Crosses: scaling of S(j c 



with S E as for multifractal G. 



S,jen, gas includes gas of local volume density at least 26po, and 
Egas gas of any local volume density. Stars: same scaling, but 
with Egas now including only gas of local volume density at least 
5po- The dashed and dotted lines are fits to the two scalings, and 
have slopes of 3.6 and 1.5 respectively. 



on the value of the dense gas density threshold, p^ en in units 
of the mean density of the cube, po- The points in the lower 
panel are generated through fitting power laws to scatter 
plots similar to those of Fig. [T] where Ed cn , ga s is, in each 
case, calculated with the appropriate definition of pden- The 
error bars in the data points represent the Icr uncertainty 
on the slope of those fits. For higher values of the threshold 
pden, fewer points are found within the cube which include 
mass at such densities, and for this reason the uncertainty 
on the scaling slope increases. From the two panels of Fig. 
[3] we can draw the following conclusions. 

(i) When paen/po is very low, most of the gas mass in 
the cube is counted as "dense gas" , and the slope n of the 
Sden,gas-S g as is close to unity. As pden increases, the scaling 
slope also increases; however, this trend does not continue 
indefinitely. As long as "dense gas" is at least one order of 
magnitude denser than the mean ISM density, the slope of 



the nonlinear scaling between Edo 



and S Eas saturates to 



the slope n (the scaling between E. 



den, gas 



and Egas, see Eq.[3j) 



some value and no longer changes appreciably with increas- 
ing Pden- 

(ii) From the lower panel of Fig. [3] we can also identify a 
trend of the saturation value of the slope to decrease as the 
highest density achieved in the multifractal cube (in units 
of the mean density) increases. The saturation value of the 
slope then settles at a value between 1.5 — 1.6 for multi- 
fractals achieving increasingly high maximum local densities 
(e.g. multifractals A, B, and C, shown with the thin solid 
black line, red dotted line, and purple long-dashed line in 

Fig.©. 

(iii) Although the maximum density achieved in the mul- 
tifractal cube is a useful first tool to discern the trends in 
the slope of the Eden, g as-E g as relation, a better quantity that 
can be used to characterize the properties of the multifractal 



6 K. Tassis 



is the width of the global PDF. As can be seen by compar- 
ing the upper and the lower panels of Fig. [3] as long as the 
global PDF is "wide" (its full width at half-maximum is 
roughly greater than 2 decs in density), the slope saturates 
to a value between 1.5 and 1.6, whereas when the PDF is 
narrower, the scaling law slope steepens and n saturates to 
increasingly high values for increasingly narrow PDFs. Mul- 
tifractal C (long-dashed purple) and multifractal D (thick 
solid green) achieve comparable maximum local densities, 
however multifractal D has a considerably narrower global 
PDF and as a result n saturates to higher values. 

(iv) Physically, a narrow PDF corresponds to a "puffy" 
ISM, characteristic of lower-mass systems with shallower 
gravitational potentials and shallower density gradients (e.g. 
Wada and Norman 2007; Tassis, Kravtsov and Gnedin 
2006). An increase of n (or, equivalently, n gas ) in low mass 
systems such as the one seen in Fig. [3] is consistent with ob- 
servations of dwarf galaxies (e.g. de Blok and Walter 2006; 
Heyer et al. 2004; Leroy et al. 2006). 

We discussed above the dependence of n on the density 
threshold pd on relevant to Eden, gas- We would now like to ad- 
dress how n changes if E gas no longer represents the total gas 
surface density, but only includes part of the gas, and more 
specifically gas of density higher than some threshold p gas - 
There are observational indications that the value of n gas 
(and thus of n) depends on p gas . For example, observations 
of dwarf galaxies indicate that when the star formation sur- 
face density is plotted against the surface density of denser 
(e.g. molecular) gas, rather than against the total gas sur- 
face density, the scaling has a smaller slope (e.g. Heyer et 
al. 2004, Leroy et al. 2005, 2006). There is a simple intuitive 
explanation for this result. As p gas increases, the densities 
of the gas included in E gas and Ed on , ga s become increasingly 
similar, and so the scaling shifts to slopes closer to linear. 
In the extreme case where p gaa becomes equal to pden, the 
scaling between the two quantities becomes identically lin- 
ear. 

We explicitly demonstrate that this trend is reproduced in 
multifractal ISM geometries in Fig. 3] We plot, for multifrac- 
tal G (a multifractal with a "narrow" PDF, appropriate for 
dwarf galaxies), Ed cn , gas against E gaa in two different ways. 
In the first case (depicted by crosses), gas of all densities 
is included in E gas (p gas is zero). The best-fit power law to 
this scaling has a slope of 3.6. In the second case, only gas 
of volume density higher than five times the mean density 
in the cube is included in E gas (p gas = 5po). As expected, 
the best-fit power-law to this scaling is shallower, and has a 
slope of 1.5. In both cases, the density threshold for inclu- 
sion in Edcn, g as is 26 times the mean density of the cube. We 
should note that the two threshold densities for E gaa are not 
special and are simply used as examples of the intuitively 
expected trend of the slope. 



4 DISCUSSION 

It is noteworthy that a geometrical interpretation for the 
Eden, gas — Egaa scaling has an interesting consequence for the 
star formation law itself. Observations of Esf suggest that 
the scaling of Esf with Ed en , g as is linear (Gao & Solomon 
2004, Wu et al. 2005), which implies that the star forma- 
tion timescale for very high density gas is roughly constant 



(see Eq. [4} and independent of Eden, gas- This however in 
turn implies that Esf scales with E gas in the same way as 

Eden, gas ■ 

Esf oc E gas . (8) 

Thus, if we assume that (a) the ISM has a multifractal topol- 
ogy and (b) the star formation timescale at high densities is 
roughly constant and does not depend on Edcn.gas, we arrive 
to a geometrical interpretation of the star formation law. It is 
tantalizing that if we do accept these two assumptions, many 
properties of the star formation law are recovered naturally. 
The robustness of the 1.5 slope for many different systems is 
the result of 1.5 being the saturation value of n for all wide 
global PDFs, which are expected for high-mass systems. The 
steepening of the star formation law in dwarf galaxies (de 
Blok and Walter 2006; Heyer et al. 2004; Leroy et al. 2006) 
is the result of a higher saturation value of n for narrower 
PDFs, appropriate for low-mass systems (Wada & Norman 
2007). The dependence of the star formation law slope on 
the tracer used to determine E gas , and the decrease of the 
slope for tracers sensitive to denser gas (e.g. Heyer et al 2004, 
Leroy et al. 2005; 2006), is the result of the dependence of n 
on the density threshold of gas included in E gas , and its de- 
crease when this threshold increases (see discussion of Fig. 

Thus, in this framework, both the universality of the 
law, as well as deviations from it under certain conditions 
and in specific systems, have a uniform interpretation as 
effects of the geometry of the ISM. Because multifractal ge- 
ometries are the ubiquitous outcome of the combined effects 
of a variety of nonlinear processes such as the ones shaping 
the scalings of the ISM, this interpretation is not based on 
any assumption about which energy inputs and dynamical 
processes dominate in ISM physics. Conversely, reproduc- 
tion of the star formation law by any ISM model cannot be 
on its own regarded as an indication that the model is a 
complete or accurate description of ISM dynamics; rather, 
it indicates that the model succeeds in reproducing a mul- 
tifractal geometry with a wide enough PDF from which the 
star formation law comes about as a natural outcome. 

The connection between a lognormal PDF and the star 
formation law has been suggested by Elmegreen (2002b), 
however in this case a complicated dependence of the effi- 
ciency on density above the threshold was also introduced, 
in contrast to our work, where no such dependence is either 
assumed or found to be necessary. Results from cosmologi- 
cal simulations (Kravtsov 2003) also point in the direction 
of the star formation law arising from lognormal PDFs with 
a constant small-scale star formation timescale, consistent 
with the picture suggested by the observations of Gao & 
Solomon (2004) and Wu et al. (2005). 

From a theoretical point of view, a roughly constant 
timescale for high-density gas can be attributed to star for- 
mation being a threshold phenomenon, in the sense that 
stars form only out of gas in the high-end tail of the PDF 
which exceeds some threshold in local volume density. Ob- 
servationally, this concept is supported by studies of the 
local environments of protostars, which are always found in 
the densest parts of molecular clouds, the dense molecular 
cloud cores (e.g. Enoch et al. 2007). Once the local volume 
density exceeds some value such that gravity overtakes the 
forces supporting the star-forming cloud, the dense gas col- 



The Star Formation Law in a Multifractal ISM 7 



lapses to form stars within a timescale corresponding to the 
threshold density. The nature of the relevant threshold and 
associated timescale depends on the relative dynamical im- 
portance of different forces in star-forming clouds. For ex- 
ample, in magneticly supported clouds, the threshold (that 
needs to be reached through ambipolar diffusion) is the crit- 
ical mass-to-flux ratio, while the relevant collapse timescale 
once the threshold is exceeded is the magnetically diluted 
dynamical timescale (see e.g. reviews by Mouschovias 1987, 
1996 and references therein). If magnetic fields are not dy- 
namically important, the relevant timescale (free-fall time) 
is very similar (see e.g. Vazquez-Semadeni et al. 2005 and 
references therein) . Although these timescales depend on lo- 
cal density, they have a fixed value tsf, dense at the threshold, 
and this may be the origin of the Gao and Solomon (2004) 
result. 

Even if the detailed local conditions in individual molec- 
ular clouds (magnetic field, Mach number distribution) differ 
so that the threshold density, as well as the associated star 
formation timescale at the threshold, differ as well, the star 
formation law will not be affected as long as the variations 
are not too large . The threshold density for star formation 
is determined by the structure of the molecular cloud cores 
(at scales of ~ 0.1 pc). On the other hand, Eden, gas is de- 
termined by the masses, numbers, and distributions of these 
cores within their large-scale environment (~ kpc), rather 
than their density structure. As a result, the star forma- 
tion threshold density and Ed C n,g as are independent quanti- 
ties. Consequently, tsf will also be independent of Egas.den.. 
Therefore, variations in tsf will increase the scatter of the 
scaling between Esf and E ga s,den but it will not change its 
slope. As long as the variations and the associated scatter 
are not so large that the scaling is lost, the star formation 
law would not be affected. 

Krumholz and Thompson (2007) have argued that the 
reason behind the Gao and Solomon (2004) findings is that 
the gas density tracer used in the particular study is sensitive 
to only a very narrow range of densities, which corresponds 
to a specific dynamical timescale value, while the star for- 
mation timescale in general does depend on the density, so 
sampling a wide density range results in the nonlinearity 
of the star formation law. Although this explanation may 
reconcile the traditional picture with observations of n gas 
between 1.0 and 1.5, it cannot account for values greater 
than 1.5, such as those suggested by observations of dwarf 
galaxies. 

Our analysis was performed, as a first step, for a three- 
dimensional cube, rather than a disc geometry which would 
be more intuitive and natural for the description of star- 
forming galaxies. As a result, the length along the line of 
sight of the cube fragments over which surface densities were 
averaged were generally larger than their plane-of-the-sky 
dimensions. On the other hand, in observations quoted here, 
the line-of-sight length of the regions studied were generally 
comparable to their plane-of-the-sky dimensions. We have 
tested that the dependence of our scalings on the plane-of- 
the-sky size of the "patches" we have used is minimal. How- 
ever, we plan to return to this problem in a more detailed 
analysis for different ISM geometries in the future. 

We should finally note that in our models, the multi- 
fractal properties are set by the mass fractions /;. In turn, 
the values of /, represent a parametrization of the shape 



of the multifractal PDF. The values of fi for the specific 
models discussed in this work were chosen so as to represent 
a variety of PDF shapes appropriate for a range of possible 
ISM geometries. The correspondence of specific PDF shapes 
to particular astrophysical systems is left to be investigated 
through simulations and observations. 



5 SUMMARY 

Observations of the scaling of Esf with the surface density 
of gas of different local volume densities indicate that a non- 
linear scaling also exists between the surface densities of very 
dense gas (traced e.g. by HCN) and total gas in star-forming 
galaxies, Eden.gas Eg l as . In this paper, we have explored the 
possibility that this nonlinear scaling may arise as a result 
of the multifractal topology of the ISM. We have investi- 
gated the properties of this scaling in simulated multifractal 
cubes, the geometrical properties of which are characterised 
by the values of the mass fractions /; that each cell inherits 
from its parent cell and the number of levels of refinement, 
H. Our findings can be summarized as follows. 

• A nonlinear scaling between the surface densities of very 
dense gas and total g as, Eden, gas oc Eg as is a natural prop- 
erty of multifractal geometries. The slope n of the scaling 
generally depends on the minimum density of the gas that 
is considered as "very dense", and tends to increase as this 
minimum density increases. However, this trend does not 
continue indefinitely, and n saturates to some value once 
the "very dense gas" threshold becomes higher than about 
one order of magnitude above the mean density. 

• The saturation value of n depends on the width of the 
global PDF. The saturation slope decreases as the width of 
the PDF increases, but again this trend is not continued 
indefinitely, and the saturation slope settles to a value be- 
tween 1.5 — 1.6 for PDFs with full-width at half-maximum 
larger than about 2 decs in density. The robustness of this 
scaling for wide multifractal PDFs is reminiscent of the uni- 
versality of the scaling slope n gas (« n, see Eq. [3J) of the 
star-formation law. 

• The origin of the nonlinearity of the Eden, gas — E gas scal- 
ing is the property of the local PDFs to extend to higher 
volume densities as the total surface density increases. This 
variation in spread between local PDFs is an intrinsic prop- 
erty of multifractal geometries. 

• With the additional assumption that, at very high den- 
sities, tsf is roughly constant and independent of Eden, gas 
(consistent with observations by Gao and Solomon 2004 
and by Wu et al. 2005), a geometrical interpretation of the 
Eden, gas — E gas can be extended to the star formation law 
itself. Since in this case Esf oc Eden, gas, we immediately ob- 
tain Esf oc Eg as . It is intriguing that in this framework two 
assumptions (multifractality of the ISM and roughly con- 
stant tsf at very high densities) are sufficient to explain 
many different features of the star formation law, includ- 
ing the robustness of the 1.5 slope across different systems 
and morphologies, the increase of the slope in dwarf galax- 
ies, and the decrease of the slope when tracers sensitive to 
higher densities (e.g. CO) are used. 
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